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BIOGRAPHICAL 


The  vibrational  modes  of  CO  adsorbed  on  a Ni(001)  metal  surface 

A short-time  expansion  is  introduced  for  the  slov  modes  of  the 

transferred  into  the  slov  modes. 

A He/OC-Ni (001 ) potential  is  constructed  from  a strongly 

molecule  and  a van  der  Uaals  attraction  accounting  for  the  He 
interaction  vilh  the  free  Ni(OOl)  surface.  A potential  is  also 
presented  for  the  Li ‘/OC-Ni (001 > vhere  a coulombic  term  is  introduced 
to  account  for  the  image  force. 

Trajectory  studies  are  performed  in  three  dimensions  to  obtain 
effective  classical  cross-sections  for  the  He/0C-Ni(001)  and  L1*/0C- 
Nl(OOl)  systems.  These  cross-sections  are  then  used  to  calculate  the 
double  differential  cross-sections  per  unit  solid  angle  and  energy 
transfer  for  the  same  systems.  Results  for  the  double  differential 
cross-sections  are  presented  as  functions  of  scattering  angles, 
energy  transfer  and  collisional  energy.  Temperature  dependence 
results  are  also  analyzed.  Extensions  of  the  approach  and  inclusion 

applications  of  the  approach  to  higher  coverages  are  discussed. 


In  the  past  few  years,  there  has  been  an  increasing  interest  in 
the  study  of  metal  surfaces  by  a wide  variety  of  experimental 

surface  reactions,  or  catalysis. 

deriving  information  about  changes  in  molecular  bonding  which  occur 
sites  (Rocca  et  ml.,  1986). 

Vibrational  spectra  of  adsorbates  have  been  observed  with 
infrared  reflection  spectroscopy  (IRS)  and  with  low-energy  electron 
scattering  (EELS).  These  techniques  are  in  practice  sensitive  only 
to  the  dipole-active  high  frequency  modes  of  vibration,  which 
correspond  to  either  intrinsic  vibrations  of  the  free  molecule,  or  to 

surface.  A further  class  of  low  frequency  modes,  the  hindered 

and  the  limited  spectral  range  and  sensitivity  of  IRS.  A technique 
which  is  in  principle  an  ideal  method  for  studying  these  vibrations 
is  inelastic  atom  scattering  spectroscopy.  It  also  has  the 


capability  to  observe  the  dynamics  of  surface  adsorbate  layers 


sensitive 


Ni(OOl)  (Ibanez  et  al.,  1983)  and  Be  vith  CO  on  Pt(lU).  In  the 
latter  system  cross-sections  have  been  measured  as  a function  of  the 
velocity  of  the  incident  He  atoms  and  the  angle  of  incidence 
(Poelsema  et  al-,  1983)  and  as  a function  of  scattering  angle  and 
momentum  transfer  (Lahee  et  al.,  1987). 


energies  in  the  thermal  regime.  Ions  have  also  been  used  as  probes 
in  experiments  of  scattering  by  adsorbates.  In  one  of  the  pioneer 
works  in  the  area  (Hulpke,  1975),  energy  and  angular  distributions  o: 
Li*  ions  scattered  by  clean  W(110)  surfaces  and  by  adsorbates,  such 
as  0-  and  CO,  were  measured  for  beam  energies  between  2 and  20  eV. 

Scattering  intensities  in  the  systems  Li‘-0/Ni(110)  and 
He*-0/Ni(110)  have  been  measured  as  a function  of  the  scattering 
angles  and  relative  energies  of  the  ions  (Bnglert  et  al.,  1983). 
Recently,  measurements  on  angle  and  energy  distributions  of  K ions 
scattered  from  V(110)  in  the  energy  range  of  12  to  100  eV  have  been 
published  (Tenner  et  al.,  1986).  T 


scattering  by  adsorbates  have  concentrated  on  elastic  scattering. 
Furthermore,  since  a detailed  theoretical  treatment  is  virtually 
impossible  to  accomplish  because  of  the  computing  time  and  memory 


required,  approximate  schemes  ate  a necessity.  Several  approaches  u 
the  scattering  of  atoms  with  thermal  energies  by  adsorbates  have  beet 


le  site  occupied  by 
Has  at  low  coverage 


based  on  extensions  of  methods  vhi 

scattering  by  overlayers  (Levi,  1982).  In  thi: 

substrate  and  does  not  contain  information  on 
the  adsorbate.  A theory  for  scattering  by  ada 

adatom  is  simply  given  by  an  effective  scattering  amplitude  obtained 

phase  scattering  amplitude  for  the  adatom.  This  theory  has  been 
applied  to  the  system  He-OC/Pt  using  the  distorted  wave  Born 

been  calculated  from  the  formalism  of  close  coupled  equations. 

computational  complexity  of  the  problem.  Quantal  approaches  have 
been  based  on  vavepackets  (Gerber  et  al.,  1984)  and  on  multiple- 
collision  expansions  (Singh  et  al.,  1986).  With  a fev  exceptions 
(Vilallonga  and  Rabitz,  1986)  theoreticians  have  considered  only 
elastic  Bragg  scattering  which  gives  information  on  the  structure  of 

scattering  from  adsorbates.  However,  the  inelastic  scattering  is  of 
great  interest  because  it  involves  the  excitation  of  surface  phonons 
and  the  vibrational  modes  of  the  adsorbed  molecule. 


scattering  by  isolated  molecules 
differential  cross-sections,  per  unit  < 
transfer.  Furthermore,  the  quantities 

distributions  o 


instead  one  vants  double 
statistical 


ates  for  the  given  surface  t( 
adsorbed  on  a metal  surface. 

Recent  experiments  have  shovn  (Berndt  et  al.,  1987)  that 


molecules  as  well  as  on  the  vibrational  modes  of  adsorbed  molecules 
can  be  obtained  from  isolated  adsorbates.  Therefore,  ve  are 
interested  in  a model  vhich  would  be  applicable  to  low  coverages,  and 
in  a simple  theory  vhich  could  provide  insights  into  the  type  of 
interactions  and  also  be  able  to  predict  trends  in  the  energy  and 

Standard  treatments  of  molecular  collisions,  based  on  target 

are  usually  unknown.  Instead,  one  can  use  the  formalism  of 

appropriate  for  scattering  by  many-body  targets  and  includes 
statistical  averages. 

strong  short  range  repulsive  forces  and  lasts  a short  time  compared 


with  the  periods  of  internal  motions  of  the  target,  so  that  an 
impulsive  model  is  appropriate.  Furthermore,  as  long  as  the 
projectile  is  light  compared  vith  target  atoms,  only  a single 
collision  occurs  vith  the  adsorbate  before  the  projectile  flies  avay. 
These  tvo  features  lead  to  a simple  expression  for  the  time- 
correlation,  vhich  can  then  be  obtained  analytically  from  the  normal 
vibrational  modes  of  the  target. 

The  vibrations  of  an  adsorbed  molecule  can  be  modelled  by  a 
cluster  vhich  Includes  the  molecular  atoms,  the  surface  atoms  binding 
the  molecule,  and  their  neatest  neighbors.  The  normal  modes  of  the 
target  can  be  grouped  into  tvo  sets,  one  vith  high  frequencies  (fast 
modes)  similar  to  those  of  the  isolated  molecule  and  the  other  vith 
lov  frequencies  (slov  modes)  characteristic  of  the  surface 
vibrations. 


When  the  target  involves  slov  modes,  one  can  approximate  the 
time-correlation  functions  by  short  time  expansions  vhich  lead  to 
Gaussian  distributions  of  the  energy  transferred  into  the  slov  modes. 

formalism  of  collision  time-correlation  functions  is  presented 
together  vith  expressions  for  the  differential  cross-sections  within 
the  impulsive  model,  and  its  application  to  scattering  by  adsorbates. 
Tvo  models  are  developed  in  Chapter  3 to  represent  on-top  and  bridge 
adsorption  sites  of  a CO  molecule  adsorbed  on  Ni(001).  Depending  on 
the  angle  of  scattering,  the  interaction  between  the  projectile  and 
the  target  may  involve  only  the  single  collision  vith  the  adsorbate, 
or  that  folloved  by  a collision  vith  the  solid  surface.  In  Chapter  4 
ve  develop  a statistical  model  to  calculate  contributions  of  double 
collisions  to  the  cross-sections.  In  order  to  calculate  the  double 


differential  cross-sections,  one  first  needs  to  obtain  effective 
the  separate  adsorbates  do  not  overlap  and  that  the  translational 

potential  for  the  systems  Li*-0C/Ni(001)  and  He-0C/Ni(001  are 
constructed  and  trajectory  studies  are  performed  to  obtain  the 
effective  classical  cross-sections.  In  Chapter  6 results  of  double 
differential  cross-sections  for  the  systems  Li *-0C/Ni (001 ) and  He- 
0C/Ni(001)  are  given  as  a function  of  energy  transfer  and  scattering 
angles.  Finally,  in  Chapter  7 the  advantages  and  disadvantages  of 
the  approach  are  presented  together  uith  some  suggestions  on  possible 
improvements  of  the  approach. 


COLLISIONAL  TIME- CORRELATION  FUNCTION  APPROACH 
TO  ENERGY  TRANSFER 


from  the  projectile  into  the  target  or  out  o£  the  target  to  the 
projectile.  The  amount  of  energy  transferred  is  determined  to  a 

it  is  expected  that  atomic  motions  within  the  target  will  place 


The  theoretical  description  of  the 
and  extended  targets  can  be  done  using 


mechanical  methods  are 


usually  based  on  expansions 

into  account  the  quantal  nai 

such  an  approach  the  target  is  described  as  a many-atom  syste 
atom-pair  correlation  functions  of  the  target  play  an  imports 
in  the  description  of  the  process.  In  this  chapter  we  review 
main  features  of  this  many  body  theory  in  its  exact  form  and 
discuss  a simplification  of  the  approach,  namely  the  impulse 
approximation,  which  makes  the  approach  more  tractable  for  il 
application  to  the  scattering  of  atom 


:he  large  number  of  energetically 
;.  An  alternative  approach  which  takes 

i been  developed  (Kioha,  1979a)  and 

1979b).  In 


hy  adsorbates. 
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, and  they  only 


) . (2n/K)  £ |<  ic*  v'  |T | kv 
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: the  completeness 
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(2.9a) 

x exp(iRxt/H)  T exp(-iHy t/K)  | lev  > , 
where  Hy  is  the  target  Hamiltonian,  and  z is  the  energy  transfer 
defined  by, 

c „ M2|k2  - (k‘)2)/2M  , (2.9b) 

The  operator  T depends  on  all  variables,  so  its  matrix  elements 
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dynamical  approximation  to  the  transition 
vay , as  it  is  shown  in  the  following  sect: 


: in  a systematic 


In  the  regime  of  hyperthermal  collisions,  the  projectile  probes 
the  internal  regions  of  the  interaction  potential  V which  is  of  a 

(Hlcha,  1979a) 


2.13) 


atom  pair  and  depends  on  the  electron  distribution  of  A and 

In  the  formalism  of  multiple  scattering  (Rodberg  and  Tib 

expressed  in  terms  of  the  two  body  potentials  v , leading  t< 
final  channel  decomposition 

T - E T(a>, 


e projectile  atom  A and  the 


Introducing  the  assumption 

calculate  T , a N * 1 body  problem  must  still  be  solved  because  G- 
involves  the  motions  of  all  the  atoms  in  the  target.  This  may  be 
substantially  simplified  whenever  the  energy  transfer  occurs  in  the 
course  of  an  impulsive  collision  of  the  projectile  A with  atom  a. 

The  basic  assumption  of  the  impulse  approximation  is  that  a large 

target  through  the  Interaction  between  the  pair 
(A, a),  while  the  remaining  N-l  atoms  in  the  target  provide  the 
restoring  forces  on  a which  determine  its  momentum  distribution 


transferred  tc 


) pair-Hamiltonian 


(2.15a> 

(2.15b) 


!.15c) 


Ka  . -h2^  /(2»a)-  , 


VE>  ’ *0,a(E)  * *0,a(E)  <HX  * Eo  - Ka>  G0(E>  • <2,17> 

Substituting  Bq.  (1.17)  in  the  equation  for  T , one  finds 

Ta  - ‘a  * 'a^O  “ *o,a>  Ta  1 (2a8) 

ta(E)  . va  . va  g0  a(E)  ta(E)  . (2.19) 

Provided  G-  and  g-  are  close  in  regions  vhere  t is  different  from 
zero,  one  can  drop  the  second  term  in  Eq.  (2.18)  and  write  Ta  = tfl(E) 
which  constitutes  the  impulse  approximation. 

R(it  -»£').  £ R(l>a)  , (2.20a) 

«<ba>  -1  S < k'V  |t„|  kv  >* 


which  shovs  the  contribution  of  each  (ab)  atom  pair  in  the  target. 
Double  differential  cross-sections  may  now  be  obtained  (Micha,  1979b) 
from  the  Eqs.  (2.20)  and  one  finds 


I(W  . 2STb  <*-,(?>*  ,(*->  (It) 


second  factor,  S<ab\iU,  is  the  Fourier  tran 


(2.21b) 
|u  > . 
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&5-  E 8>  S1"1  Ao  ■ «.!*) 


multiple  collision  terms  would  be  small  for  mA/mx  « 1 and  kR^  » 
interatomic  distances  in  the  target.  Geometric  considerations  will 
general  the  validity  of  these  assumptions  mus 


I Application  to  Scattering  by  Adsorbates 


when  adsorbates  are  introduced, 
the  loss  of  periodicity  which  i 


crystal  surface  can  occur 
associated  with  bare  crystal 


ie  adsorbates  form  a regular  lattice  01 


cell,  usually  larger  tha 


e original  ui 


may  think  of  a single  adsorbate  or 
arrangement  of  adsorbates.  The  fi 

the  periodicity  of  the  system  as  t 
surface,  is  lost,  usually  because 
the  atom  spacing  in  the  substrate. 


the  surface  is  destroyed.  Here  we 
>n  the  surface,  or  a non-periodic 


tie,  parallel  to  the  crystal 


adsorbates  requires  a variety  of  approaches,  depending  on  the 


magnitudes  of  collision  energies  and  on  the  extent  and  type  of 

The  many-body  approach  has  been  previously  used  to  study  energy 
transfer  in  hyper  thermal  collisions  of  ions  with  a solid  surface 
(Micha,  1981).  In  such  cases  the  scattering  intensity  is 

single-atom  collisions  contributing  over  all  scattering  angles  9 

For  scattering  by  light  projectiles  with  energies  in  the 

of  the  adsorbate  and  of  the  atoms  in  the  uppermost  layer  of  the 

similar  way  as  in  the  case  of  surface  scattering 


ection  for  the  projectile-a 
ie  on  the  plane  (x,y,z»0)  a: 


:et,  which  can  now  be  adsorbate  atoms 
s the  elastic  differential  cross- 

ie  adsorbate  is  placed  at  a distance 


Z , given  by  the  equilibrium  position  of  the  adsorbate  above  the 
surface  plane.  The  projectile  moves  in  the  direction  of  decreasing 


state  of  atom  1 in  the  target  and  can  be  obtained  from  semiempirical 
potentials  whose  parameters  have  been  fitted  to  reproduce 
experimental  data. 


elastic  differential 


S<11) 


has  the  meaning  ol 


energy  transfer  e. 

This  self-correlation  function,  must  be  obtained  from  a 

model  Hamiltonian  for  the  system.  A procedure  which  can  be  used, 
appropriate  in  cases  of  lov  coverage,  is  to  describe  the  target  by  a 

the  vibrational  modes  which  are  localised  in  the  vicinity  of  the 

surfaces,  it  is  expected  that  double  collisions  may  contribute  at 
certain  scattering  angles.  The  application  of  Eq.  (2.25)  is  only 
valid  in  regions  of  single  collisions.  However  a statistical 
approach,  which  uses  the  geometry  of  the  problem  and  allows  for  the 
description  of  double  collisions,  will  be  presented  in  a later 
chapter. 


CHAPTER  3 

VIBRATIONAL  FREQUENCIES  AND  NORMAL  MODE! 
OP  MOLECULES  ADSORBED  ON  SURFACES 


ie  studied  with  a variety  of  experimental  techniques  involving 


scattering  of 

adsorbate  and 
determines  the 


electrons  (Ibach  and  Mil 

its  substrate,  and  about 
probing  particles  are  ate 


, 1982),  photons,  neutrons 
out  the  conformation  of  the 


d reshaping  of  spectral 


',  possibly  resulting  ii 

i be  simplified  when  the  internal 
identified  and  separated  from  the 


interpretatioi 

The  normal  modes  of  adsorbates  on  solid  surfaces  can  be  obtained 
from  generalizations  of  theories  developed  for  surface  lattice 
dynamics.  Approaches  can  be  based  on  Green  functions  (Rahman,  Black 
and  Mills,  1982),  slab  models  vith  a small  number  of  layers  (Strong 
el  al.,  1982),  and  cluster  models.  Green  functions  require* exte 
work  vhile  slab  calculations,  although  less  involved,  may  not  a)1 
properly  describe  the  influence  of  site  symmetry  (Rahman  et  al.. 
1983).  Furthermore,  they  have  been  used  together  vith  simple  foi 
fields  of  the  central-field  type.  However,  one  expects  that  mnn 
detailed  force  fields  will  be  needed,  containing  In  particular 


excitation  of  adsorbate  vibrations. 

Cluster  models  have  been  found  adequate  to  describe  adatoms 
(Black  et  al.,  1982),  and  in  some  cases  have  been  as  effective  as 
models  based  on  periodic  small  clusters  (Black,  1982).  Applied  to 

vhat  is  knovn  about  the  vibrational  modes  of  the  isolated  molecule 
adsorption  site. 

Experimental  studies  (Andersson,  1977)  of  CO  on  Ni (001 ) have 
provided  evidence  that  CO  molecules  can  bind  on  a cop  site  above  a Ni 
atom,  on  a bridge  site  betveen  tvo  of  them  or  in  a mixture  of  both 
(Bertolini  and  Tardy,  1981).  For  a dilute  CO  lattice  gas  on  Ni(001), 
tvo  peaks  are  observed,  at  239  and  256  meV,  in  electron  energy  loss 
spectra  (EELS).  Hovever,  at  higher  coverages  only  a single  peak  at 
256  meV  is  observed.  These  findings  have  been  interpreted  as  due  to 
the  presence  of  both  top  and  bridge  species  in  the  case  of  low 
coverages,  but  of  only  the  on-top  species  in  the  higher  coverages 
(Andersson  and  Pendry,  1979). 

Model  calculations  on  small  clusters  have  been  performed  by 

adsorbed  on  a nickel  surface.  The  frequencies  of  a NijCO  cluster  for 
the  on-top  site  and  a Ni -CO  cluster  for  the  bridge  site  have  been 
analyzed  (Richardson  and  Bradshaw,  1979),  in  terms  of  semi-empirical 
force  fields,  with'  only  the  C and  0 atoms  moving.  The  frequencies  of 


(Allison  and  Coddard,  19 

adsorption  and  apply  it  to  CO  adsorbed  on  Ni(OOl).  It  is  based  on  a 

literature,  and  it  incorporates  site  symmetry.  It  also  allows  for 
moving  surface  atoms  so  that  a more  detailed  picture  of  the  adsorbate 
dynamics  can  be  derived.  Furthermore,  the  substrate  cluster  contains 
two  surface  modes  whose  frequencies  can  be  fitted  to  known  values  at 

1955).  The  numerical  results  and  figures  shoving  the  normal  modes 
are  given  in  section  3.3. 


3.1  Cluster  Hodels  and  Force  Fields 
To  model  adsorption  at  the  on-top  site,  we  use  an  eleven-atom 
cluster,  shown  in  Fig.  3.1,  containing  C,  0 and  nine  Ni  atoms,  with 

metal  atoms  that  are  nearest  neighbors  to  the  bonding  site,  four  in 
the  first  layer  and  four  in  the  second  layer,  with  the  CO  bonded  to  a 

In  the  bridge-bonded  situation,  we  have  chosen  a 16-atom  cluster, 
shown  in  Fig.  3.2,  that  contains  fourteen  Ni  atoms,  of  which  the  two 
bonded  to  C are  alloved  to  move.  This  cluster  also  contains  all 


correspond  to  surface  atoms,  and  vith  numbers  6 to  9 t 
nearest  neighbors  in  the  second  layer;  C,0  and  nickel 


Figure  3.2:  Bridge  Cluster.  Nickel  atoms  with  numbers  1 to  8 

correspond  to  surface  atoms,  and  with  numbers  9 to  Id 
nearest  neighbors  in  the  second  layer;  C.O  and  nickel 

height  of  the  C atom  above  the  surface  (d„)  is 
determined  by  the  angle  Ni-C-Ni. 


representing  the  bonding  site,  six  in  the  first  layer 
second  layer.  The  carbon  atom  of  the  CO  molecule  is  fc 
nickel  atoms,  one  in  a position  corresponding  to  a fac 


A C^v  point  group  can  be  associated  to  the  top  bonding  situation 

Since  on  a real  surface  all  translation  and  rotations  are  restricted, 
we  expect  all  3N  degrees  of  freedom  to  be  vibrational  modes. 
Therefore,  we  should  find  3 x 3 « 9 vibrational  modes  for  the  top 


be  classified 


s and  four  modes  which  are  equivalent  to  the 
of  an  XYZ  isolated  linear  molecule  in  the  gas 


The  latter  could  also  be  considered  as  arising  from  the  vibration 
substrate  represented  by  a single  nickel  atom. 

of  these  correspond  to  the  molecular  modes  of  a X-CO  isolated 
molecule  in  the  gas  phase.  Again  this  could  also  be  considered  as 
the  vibration  of  the  CO  molecule  interacting  vith  four  phonon  modes 
of  the  substrate  represented  by  two  nickel  atoms. 


the  actual  form  of  the  normal 
vibrational  modes  according  tt 


the  symmetry 


classify  the  different 
he  potential  energy  of  the 


clusters. 


Usually,  the  potential  energy  is 


internuclear  distances,  one  each  for 
distances,  one  for  the  change  in  the 


and  in  a number  of  angles.  For 
f 12  internal  coordinates,  eight 
changes  in  the  Ni-Ni 


adjacent  Ni-Ni  bond. 

Denoting  by  s.j  the  changes  in  the  Ni-Ni  internuclear  distances, 


and  s^  the  changes  in  the  angles  described  abo 


In  the  above  expression  the  subscript  1 had  been  used  to  denote  the 
metal  atom  bonded  to  the  CO  molecule,  f^P  and  f^j  denote  the  force 

layer  and  that  between  first  and  second  layer  atoms. 


it  changes  in  i 


i internuclear  distances,  3 to 


(3.2) 


the  Ni(OOl)  surface  at  different  points  of 
sone.  Two  surface  nodes  with  frequencies  of  104.7  and  125.5  era-1 

and  f ;<*>  in  2,.  (3.2). 

from  data  for  Ni(CO),,  (Jones.  1960).  For  the  bridging  species  the 

for  a molecule  of  the  type  i^CO  (Herzberg,  1965).  which  is  the 

the  Hi  bonded  to  the  CO  molecule  have  been  taken  avay.  The  chosen 
force  constants  and  geometric  parameters  are  presented  in  Table  3.1 
and  3.2  respectively,  where  . is  the  distance  between  nearest 

the  C atom  to  the  surface  in  Fig.  3.2.  As  a preliminary  check  of  the 
model  and  the  chosen  field,  we  calculated  the  frequencies  of  an 

position  cluster.  The  results  are  within  1 cm'1,  as  can  be  seen  in 
Table  3.3,  from  those  obtained  by  an  ab  initio  calculation  on  the 
same  molecule  (Allison  and  Goddard,  1982). 

For  the  cluster  representing  the  bridge  bonded  situation,  the 
limiting  case  was  11,00,  obtained  by  replacing  in  our  cluster  the 


Internal  Coordinate 


Force  Constants 


£cq  (mdyne  A-*) 


£|J{1)  (mdyne  A"1) 
(mdyne  A"*) 
f (mdyne  A rad-2) 
f.  (mdyne  A rad"2) 
t,  (mdyne  A rad'2) 


Table  3.2  Geometric  Parameters  of  Clusters 


Geometric  Parameters  Top 


dc0  (A) 

dCNi  (A) 

dN1  (A) 
dc  (A> 


angle  (deg) 


(in  cn-') 


displacements  X.  The  transformation 


2V(s)  = 


The  relations  betveen  the  different  coordinates  are  listed  in  the 

for  the  on-top  cluster  and  Table  3.5  for  the  bridge. 

Decius  and  Cross*  1955).  The  G matrix  is  defined  by  the  kinetic 
energy  T.  Starting  with 


2T  . <$)'  O’1  | , 

where  g-1  = (U-1)1  H g_1  and  M is  the  diagonal  matrix  of  mi 
From  these  F and  G matrices,  one  obtains  the  frequencies 

<C_F)  Ak  - ^ Ak  , 


V - <2ncV2  . 


r Bridge  Cluster 


Fjj  - 2(a3 


F12  = 2a4(a3  - a4)  fa  - £C(J 


F13  . -(2/J2)  a3(a3  - a4)  £a  - (2/j2)(dc/dc(a)2  £c(|,  - (4-12)  b2fg 
Pj4  - -(2/JI)  bjdcfcm  - (2/42)  a2  (a3  - a4)  £a  * (4/42)  b2a2£ 


FJ4  . -(2/41)  a2a4fa 

F33  = a3fa  * (dC/dCNl)2  £CNi  ‘ 2<dC/dCNidM)2  £b  * 2blfg 
. (40/21)  f£2) 

- ZbjBjf  - (4/7)(l/42)  fj2) 


* (17/5)  fJJ*  . (8/7)  f*2) 


B2  Symmetry 

F11  ■ a3  fa  * <dc/dcm)Z  £CNi  *■  <40/21>  fNi’  * 2(dc/dMidCNi)2  £b 
F12  . (2/42)  a3(a2  - a3)  - (l/42)(dNidc/d2m ) £c(,. 

Fj3  » (2/42)  a3a3£a 

Fw  . (dNidc/2d2Nj)  £C[|1  * (44177)  £*2)-  a2a3£a  . (d^d,,^2^)  f b 
c)2  £c[|1  * 2(a2  - a3)2  fa 


<l/42)(dN./dc 


P24  “ " (2/^2)  a2  <a2  " *i>  fa  - ( 1/2-12)  (<,Ni/dcNi>2  £cm 
F34  - - (2/J2)  aja2fa 

F44  . a j£a  * (dNj/2dCNi)2  £CNi  * £|>/<2dCNl)  £b  * (7/5)  41* 

* (8/7)  f<2> 

A-  Symmetry 

Pn  = (8/5)  f*^  * (20/21)  f£2)  * (2/d2Ni)  fb 
B,  Symmetry 

<4/dCNi)  £b  * '<dC  * dC0>/<dC  dC0>'2  £d 

- fd  (dc  . dco)/(dc  d2„) 

- 4£b/(45  d2N.)  - £d(dc  . d£0)/(j2  dc()d2> 

£d/dC0  F23  " £d/<^2  dCdC0 

(8/5)  f£P  * (20/21)  f£2)  * £d/(2djb  * 2£b/d^N- 

The  matrices  GF-  constructed  in  this  vay  reflect  the  symmetries  of  the 
clusters,  thus  simplifying  the  calculation  and  the  classification  of 
the  modes  according  to  the  symmetry  of  the  clusters.  In  the  on-top 
case  the  resulting  GF  matrix  is  block-diagonalized  into  one  3x3 
matrix  belonging  to  the  Aj  irreducible  representation,  and  one  6x6 
matrix  containing  tvo  coupled  modes  belonging  to  the  E irreducible 


F33  ‘ 


o two  4 x 4 matrices,  one  belonging  ti 
2 irreducible  representations,  one  3 : 
d one  lxl  matrix  belonging  to  the  A 


irreducible  representation. 


3.3  Numerical  Results  for  CO  on  Ni(001) 

The  vibrational  inodes  for  the  on-top  and  bridge  clusters  are 
shovn  in  Figures  3.3  and  3.4,  respectively.  The  resulting 

coordinates  are  presented  for  the  on-top  and  bridge  clusters,  in 
Appendix  C.  We  separately  discuss  the  on-top  and  bridge  cases. 
3.3.1  The  On-Top  Cluster 

The  calculated  frequencies  along  vith  previous  experimental  ai 
theoretical  results  are  presented  in  Table  3.6  for  the  on-top 

to  the  CO  stretching  for  the  free  molecule,  2143.2  cm  * , and  two 
frequencies  v*  and  vith  values  similar  to  the  experimental 
frequencies  (of  381  cm'1  for  CNi  stretch,  and  of  461  cm  1 for  NiC 
bend)  found  in  carbonyl  compounds,  Ni(CO),  (Crawford  and  Horvitz, 


1948).  We  also  find  three  frequencies  \>3,  Uj  and  v6  vith  values 
lover  than  the  largest  frequency  for  the  bare  substrate  phonon  band, 
295  cm'1  (Black  et  al.,  1981). 

The  mode  vith  the  highest  frequency,  Vj,  is  related  to  the  free 
molecule  CO  stretch,  v2  relates  to  the  Ni  C stretching  mode,  and  the 

to  the  surface. 


0 

c 

Ni 


Figure  3.3:  Normal  Vibrational  Displacements  for  the  On-Top  Cluster. 

The  u-symbol  label  different  irreducible  representations 
for  each  symmetry  group. 
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Figure  3.4:  Normal  Vibrational  Displacements  for  the  Bridge  Cluster. 
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Figure  3.4 — continued 
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periodic  (Ni)~C0  cluster  using  green  function  methods  (Black,  1982). 
Ve  also  obtained  a value  of  102  cm-1  for  the  mode  which  involves 

higher  than  the  one  obtained  by  Black  for  the  overlayer  and  for  the 
isolated  (Ni)eCO  cluster,  could  be  associated  to  the  fairly  broad 
peak  centered  at  148  cm'1  obtained  in  the  spectral  density 
calculations  for  the  full  crystal. 

3.3.2  The  Bridge  Cluster 

Asymmetry,  three  with  Bj  symmetry,  four  with  Bj  symmetry  and  one 

Vg  , have  values  similar  to  molecular  mode  frequencies,  and  the 
remaining  v^1 , v^1,  v^1 , Vjq,  v^,  and  have  values  lover  than 

the  largest  bare  substrate  phonon  frequency  (295  cm'1).  These 
calculated  frequencies  are  presented  in  Table  3.8  along  with  results 
from  other  calculations  and  previous  experimental  results. 

The  Aj  mode  with  the  largest  frequency,  Vj,  can  be  easily 

other  two,  and  u^,  can  be  considered  as  frustrated  translations  in 
a plane  normal  to  the  surface,  with  the  last  two  modes  involving  C- 
Ni-Ni  bending.  Three  bending  modes  belonging  to  the  Bj 
representation  involve  motions  out  of  the  plane  formed  by  the  two 
nickel  atoms  and  the  carbon;  one  of  these,  v^,  can  be  considered  as  a 
frustrated  translation  and  the  other  two,  and  v^,  as  OC ( N i ) ^ 
bendings.  A summary  of  these  assignments  is  shown  in  Table  3.9. 


Table  3.8  Comparison  ofResulJs,  Bridge  Cluster, 


Experiment  Our  (Richardson  & 

(Andersson,  1977)  Results  Bradshav,  1979) 


Frequency 


” KKXh'KS'S;  BSrAE1 


"s‘sa* 


frequencies  presented  here  describe  new  nodes  for  the  bridge  cluster, 
which  have  not  been  presented  in  previous  works. 

A fairly  clear  distinction,  see  energy  diagran  shown  in  Figure 
3.5,  can  be  made  between  nodes  with  relative  high  frequencies 

fron  the  interaction  of  the  adsorbate  with  the  nodes  of  the  bare 
substrate.  This  separation  of  nodes  is  nore  evident  for  the  bridge 


and  their  frequencies  and  normal  modes  can 
in  the  study  of  atoms  scattered  by 


adsorbates. 


FREE  MOLECULE 


CHAPTER 


0 MOLECULAR  ADSORBATES 


translational  degrees  of  freedom  can  frequently  be  described 
classically  while  the  remaining  ones,  only  vibrational  degrees  of 
freedom  in  our  case,  must  be  treated  as  quantized. 

It  was  pointed  out  in  Chapter  2,  that  by  relating  inelastic 

scattering  by  extended  targets  vithout  having  to  expand  wavefunctions 
in  the  many  excited  states  of  the  target.  Inelastic  scattering 
cross-sections  for  collisions  in  which  the  relative  momentum  changes 

For  experimental  conditions  of  high  relative  ei 
projectiles,  the  interaction  is  dominated  by  strong  short  range 
repulsive  forces,  and  lasts  a short  time  compared  with  the  periods  of 
internal  motions  of  the  target,  so  that  an  impulsive  model  is 

double  differential  cross-section  is  obtained  from  TCFs  involving  the 
position  operators  of  the  target  atoms.  Ue  show  in  Sec.  4.1  how  to 
obtain  the  vibrational  TCFs  analytically  from  the  normal  vibrational 
modes  of  the  target.  Furthermore,  since  the  normal  modes  of  the 


id  light 


F<aa,{*,t>  - « exp[-l?.?a« 
« ...  » - E “v<  V|  ...  I'M 


3.3.3. 


F<aa>(it,t) 


,<0>l) 


FvjCt)  = « expl-iic  u (0)J  expliK-u (t)|  » 


‘ E hi  s 


s 0j(O  - 0j(0)  copy,  which  are  r, 


*.)  • E •«  •<  «■•  » 


FviJ(t)  = ]~[  « exp[-iit-6ajO3(0)I  x exp|iK-6aJ0j(t)]  »„.  (4.7) 
j-1 

61),  through 
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o Eq.  (4.7),  we  can  « 


'jw  ■ n i“<»  . 
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total  vibrational  correlation  function  Fy(c  in  Eq.  (4.10)  can  be 
rearranged  as  a sun  of  products  of  exponentials 

PV(C(0  * E e*P(-*nBvO  P„(2)  . (4.12a) 


Where  «v  and  n are  colunn  matrices  with  elements  Wj  and  Integers  n^, 
the  vibrational  TCF  vould  consist  of  a group  of  S-function  peaks 


4.2  Short  Time  Expansion 

The  calculation  of  the  vibrational  TCF  given  in  Eq.  (4.10)  can  be 
simplified  whenever  the  vibrational  modes  of  the  target  satisfy 
H«.  < kgT.  Therefore,  one  can  approximate  by  a short  time 

In  F(t)  . In  F(0)  * t P(0)/F(0) 

* \ t2|F(0)/F(0)  - f(0)2/F(0)2l  , (4.13) 

from  Eq.  (4.10),  the  expression  for  each  of  the  lov  frequency  modes  j 
is  given  by 


r2/*2  , 


quantities  as  given  by  Eq.  <4.11). 


Collecting  all  the  lov  frequency 
final  expression 


contributions  one  obtains  the 


In  Fv(t>  - it  £ A./H  - t2  £ r?/H2  . (4.15) 

3 3 

Caussian  distribution  peaking  at  the  energy  transfer 
c = £ A,  with  a width  F = £ F. . Both  quantities  would  increase 

3 3 

with  momentum  transfer. 


Statistical  Hodel  f 


e Collisions 


t likely 


event  is  excitation  of  both  adsorbate  and  subst 
although  deexcitation  is  also  possible  from  the  thermally  populated 
excited  states  of  the  target. 

vectors.  For  an  atom  approaching  along  the  perpendicular  to  the 
in  Figure  4.1,  vhich  has  been  drawn  using  covalent  atomic  radii  plus 


parameter  is  large,  the  first  deflection  is  smaller  than  n/2  and  the 


Trajectories  shoving  regions  of  single  and  double 
collisions.  Model  dravn  using  covalent  atomic  radii 


53 


atom  next  collides  with  the  surface  before  flying  avay  at  an  angle 
larger  than  lt/2.  Hence,  the  deflection  function  of  such  a collision 
vould  shov  two  branches  for  3 > n/2.  It  is  then  obvious  that  an  atom 
which  had  been  deflected  by  either  mechanism  could  be  detected  at  the 
same  final  scattering  angle.  Therefore,  to  obtain  the  double 
differential  cross-section  one  needs  to  add  all  contributions 
involving  the  same  amount  of  energy  transfer  and  having  the  same 
scattering  angles  regardless  of  the  mechanism  by  which  this  occurred. 
In  Chapter  2 an  expression  for  the  double  differential  cross-section, 
for  single  collisions  and  under  the  impulse  approximation,  was  shown 
to  be  given  in  terms  of  the  TCF.  Using  the  vibrational  TCF  obtained 


effective  classical  cross-section  and  the  second  has  the  meaning  of  a 
probability  per  unit  energy.  Then  equation  (4.16)  can  be  rewritten 


8 - JSSr  E V*>  • E • <*'>»> 

In  particular,  for  a scattering  characterised  by  the  energy 
transfer  e = -nhu,  we  find  the  partial  differential  cross-section  to 


S-HM 

i-|HM 


& ■ tt.in 


cluster  and  nm  for  the  metal  surface  away  from  the  adsorbate.  The 
latter  contains  the  parallel  and  perpendicular  vibrational  modes  of 

(dP/dc)(II)  . | e'ie,/M  E e Pn(e)  , (4.20a) 

Pn(E)  = J dy  P c(ye)  P m |(l-y)e|  , <4. 20b) 


vhere  En  « -£  n^KWj  and  P|;  is  the  probability  of  transferring  n 
quanta  when  the  total  energy  transfer  is  £ = EC  * £m.  The 


probability  P has  been  vritten  as  the  convolution  of  the 
probabilities  P and  P for  excitation  of  the  cluster  and  the  clean 


surface,  corresponding  to  a sequence 
integration  variable  y in  Eq.  (4.20) 


defined  by  y = EC 


(t,y) 


f”  <t,l-y> 


(4.21a) 


Fc  (t,y)  = £ exp(iE  ct/K)  P e (ys)  , (4.21b) 

F"(t,l-y)  ■ E exp(iE  mt/W)  P Bl(l-y)El  . (4.21c) 


Here,  Fc(t,y)  and  Fm(t,l-y)  are  the  vibrational  TCFs  for  the  cluster 


and  clean  metal  respectively.  Each  of 
calculated  from  the  expressions  in  Eq. 


Hovever,  it  is  usually  found,  as  vas  the  case  in  Chapter  3 £o 
CO  molecule  adsorbed  on  a Nickel  surface,  that  frequencies  can  be 
separated  into  high  and  lov  ones.  Thus,  for  each  lov  frequency  m< 

, 2 .(II)  d2o 


(4.22b) 


Equation  (4.22)  m; 
cross-sections  ( 
amount  of  energy  e 


CHAPTER  5 

GAS-ADSORBATE  INTERACTIONS  POTENTIALS: 
EFFECTIVE  DIFFERENTIAL  CROSS-SECTIONS 


In  previous  chapters  ve  have  shovn  that  double  differential 
cross-sections  for  atom-adsorbate  collisions  can  be  obtained  from  a 
product  of  deflection  probability  and  a target  absorption 

differential  cross-section.  In  order  to  obtain  this  effective 
differential  cross-section  one  must  first  specify  an  interaction 

interaction  potential  and  the  usual  procedure  one  can  follov  to 
obtain  analytic  representations  which  are  desirable  for  scattering 


behavior  of  the  potential,  or 
potentials. 

helium  atoms  and  lithium  ion: 

These  model  potentials  a: 
two  systems,  He  - adsorbed  C( 


n quantitatively  about 


the 


molecule,  t 


n the  desired  effective  differential  c: 


sections. 


The  general  features  of  gas-adsorbate  potentials  are  expected  to 
be  similar  to  the  qualitative  picture  of  the  gas-surface  potentials. 

potential  has  a strong  repulsive  part  in  front  of  the  solid.  The 

the  electrons  of  the  gas  atom  and  of  the  electrons  at  the  surface  of 
the  solid.  This  gives  rise  to  a periodic  modulation  of  the  repulsive 
part  of  the  potential  parallel  to  the  surface.  At  longer  distances 


le  attractive  pt 


There  have  been  only  a relatively  small  number  of  attempts  to 
calculate  the  physical  gas-surface  interaction  potential  from  first 
principles.  Calculations  of  He  on  metals  have  been  done  (Zaremba  and 
Kohn,  1977)  to  calculate  the  attractive  van  der  Vaals  interaction  of 
large  distances  and  the  repulsive  interaction  at  small  distances  and 


procedure  the  position  of  the  reference  plane  Z = Z-,  from  which  the 
distance  of  the  gas  atom  should  be  measured,  appears  as  an  essential 
parameter. 

More  recently,  methods  for  generating  surface  potentials  from 
surface  electron  densities  have  been  presented.  Some  authors 
(Esbjerg  and  Norskov,  1980)  have  suggested  a theory  which  gives  a 


the  surface-electron  charge  density.  However,  when  experimental 
results  are  compared  with  theoretical  predictions,  they  strongly 
disagree.  In  particular,  for  the  Ni < 110)  surface  the  corrugation 

experimental  value  (Annelt  and  Haydock,  1986) . Similarly,  several  a' 
initio  calculations  of  potentials  for  weakly  corrugated  surfaces  hav, 
not  been  very  successful,  predicting  approximately  two  times  too 


methods,  allowing  tl 


n (Batra  et  al.,  1985).  However,  semi-empirical 
e adjustment  ol 
e scattering  d; 

An  approximation  frequently  used  t 
physisorption  potential  energies  o 
the  additivity  assumption.  The  thermally  averaged  solid  surface  i 

Fourier  expansion  as  follovs: 


985).  However, 
as-surface  systems  is  bt 


v(r)  = £ vg(Z)  exp(i34s)  . (5.1) 

a 

where  r = (x,y,z)  > (iL,Z),  with  the  Z-axis  extending  in  the 

the  reciprocal  lattice  vectors  of  the  surface.  The  first  term  in  the 
expansion,  VjJ(Z),  is  considered  the  over -all  average  potential  and  is 
the  large  potential  in  perturbation  theory  which  causes  specular 
diffraction.  The  v^(Z)  for  (•  * 0 are  treated  as  perturbations  which 
cause  non-specular  diffraction.  Diffraction  experiments  give 
information  on  both  the  dimensions  of  the  uni t cell  via  the  angular 
location  of  Bragg  peaks  and  the  distribution  of  the  scattering 


form,  the  potentials  vg(Z),  including  6 = 6,  must  be  simple. 

A convenient  empirical  representation  of  the  vg(Z)  is  the 

which  represents  vg(Z)  by  a Morse  potential  and  vg(Z)  by  an 
exponential  repulsion  for  6*0. 

vg(Z)  ■ D(exp(-2«Z)  - 2exp(-«Z)|  , (5.2a) 

vg(Z)/D  > Kg  exp(-2mZ)  3*0.  (5.2b) 

vhere  only  the  four  smallest  nonzero  reciprocal  lattice  vectors  are 
used,  the  sum  over  3*0  reduces  to  four  terms  and  the  Kg  reduces  to  a 


single  value. 

The  qualitative  picture  of  the  gas-surface  interaction  we  have 
described  has  to  be  modified  in  the  case  of  ion-surface  interactions. 


the  surface  interaction  c 


mes  from  the  overlap  of  the  electronic 


the  surface  the  interaction  is  attractive  and  is  well  described  by 
the  macroscopic  image  potential,  V.ffl  = -e  /4Z.  The  image  picture, 
however,  loses  validity  when  the  particle  approaches  to  vithin  a few 


Vaals  interaction,  a parameter  may  be  introduced  to  define  the 
position  of  the  image  plane  which  improves  the  validity  of  the  image 
potential. 

Recently,  a way  to  construct  the  interaction  potential  between  an 
ion  and  a surface  of  a monatomic  solid  has  been  presented  (Mann  et 
al.,  1987).  The  K*  - V(10t))  potential  constructed  this  way  was  found 


to  agree  with  on  empirical  potential  used  to  fit  scattering  data  for 
the  same  system  (Hulpke  and  Mann,  1985). 


The  interaction  of  an  atom  with  an  adsorbate  will  be  comprised  of 
direct  potentials  (the  attractive  van  der  Uaals  and  the  repulsive 


of  particular  importance  in  the  case  of  long  range  interactions  and 

The  induced  polarization  interaction  of  closed  shell  atoms  with 
anisotropic  adsorbed  molecules  has  been  investigated  and  a theory  of 
the  long  range  van  der  Vaals  potentials  acting  between  a He  atom  and 

Gumhalter,  1987),  and  the  full  potential  has  been  given  as  the  sum  of 
the  long  range  interaction  plus  the  short  range  repulsion  for  the 

A simple  approach  to  the  full  description  of  an  atom-adsorbate 
interaction  consists  in  superimposing  an  ab  initio  atom-gas  phase 
molecule  interaction  potential  on  an  empirical  atom-surface 
potential.  This  procedure  has  been  used  in  the  system  He/OC-Pt(lll) 
(Lahee  et  al.,  1987)  and  has  been  found  to  give  good  agreement  with 

5.1.1  He/OC-Ni (001 ) Interaction  Potential 


used  in  trajectory  studies 
m-adsorbate  potential  can  b 


we  follow  the  assumption 
given  by  superimposing  an 


om-molecule  potential  oi 


»B(r)  . J de  p(s.t)  g(c)  , 


V"-"  ■ vo  *""  {■  * ir  0i]«..»} , <=-*•> 


4, 4f  -j*- 


.y.*>  ■ Bs  D exp|-6(z-z0)]h(x,y)  , 
(Z)  - -Cvv/(z-zvv)3  . 
represents  the  distance  betveen  tl 


Ra  = |R%1 

and  R is  the  distance  betwei 
CO  molecule.  In  the  case  of 
experimentally  (Andersson  am 


CO  on  Ni(001).  Rm  has  been  determined 
Pendry,  1980).  The  geometry  of  the 


plane  and  the  z coordinate  is  perpendicular  to  the  surface  and  passes 
through  the  center  of  the  adsorbed  molecule. 

Prom  electronic  calculations  of  CO  on  a cluster  of  surface  atoms 
there  is  overwhelming  evidence  that  there  is  a chemical  bond 
formation  between  the  substrate  valence  band  orbitals  and  the 
adsorbed  CO  valence  molecular  orbitals.  The  CO  5c  molecular  orbital 

metal  charge  may  be  back-donated  into  the  formerly  unoccupied  CO  2n 


geometry  used  for 
> On-top  model;  b) 
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metal-CO  interaction  (Bauschlicher,  1986).  Tills  back  donation  into 
the  2n  orbital  is  then  responsible  for  the  weakening  of  the 
molecular  bond  and  an  increase  of  the  CO  bond  length.  The  C-0  bond 
length  in  a CO  molecule  is  1.12B  A,  vhereas  in  a CO.  molecule,  it  is 
1.159  A which  is  closer  to  the  value  found  in  carbonyl  compounds  and 


In  view  of  the  above  arguments, 
parameters  given  by  the  fitting  of 
for  He  - C02  (Clary,  1982),  for  the 
the  potential.  Diffraction  data 
Ni(001)  system,  so  we  have  used  the 
for  the  He  - surface  Interaction, 
have  used  in  this  work  is  given  in 


e have  chosen  to  use  the 

SCF  interaction  potential  energy 

e - adsorbate  repulsive  part  of 

arris  and  Liebsich  parameters 


considers  the  bridge  bonding  situation  is  tl 
of  the  clusters  presented  in  Chapter  3. 


which  changes  when  one 


interactions 


distances,  because  of  electronic  overlap  of  wave  functions.  In  order 
to  accomplish  this  damping,  we  have  used  a function  of  the  type  f(z) 

- exp| S(z-Zq) ]/(l  * exp | 6(2-2.) I ] which  for  z<z»  goes  smoothly  to 


The  interaction  potentials  for  the  on-top  an 
shovn  in  figures  5.2  and  5.3,  respectively.  It 


the  region  of 
are  interested.  The 


Interaction  Potential 


He/0C-Ni(001) 


A3/(Z'BH)A 


the  bridge. 


shorter  distance 


the  CO  molecule  to  the  surface  in  the 
seen  in  the  inset  of  both  figures  that 
a deeper  veil  than  the  top  potential, 
part  of  the  interaction  betveen  the  He 
stronger  when  the  He  atom  is  over  the 


This  is  because  the  attractive 


bridge  model,  the  adsorbate  is  positioned  in  betveen  tvo  surface 


(y,z)  plane  is  shovn  for  the  tvo  models.  Again  tl 
larger  turning  points  than  the  top  model  but  nov 


potential  in  the 
le  bridge  gives 


le  fact  t 


5.1.2  Li  /OC-Ni(OOl)  Interaction  Potential 

The  ion-adsorbate  interaction  potential  can  in  general  be  vritten 

interaction,  the  ion-surface  interaction,  and  the  image  interaction. 

interaction  potential  betveen  an  ion  and  the  surface  of  a monatomic 
solid  (Hann  et  al.,  1987)  is  found  to  agree  with  an  empirical 
potential  used  to  fit  scattering  data  for  the  system  Li*/tf(001) 
(Hulpke  and  Hann,  1985).  Me  have  chosen  to  represent  the  Li+-Ni  (001) 
potential  by  a modified  MEER  potential  vhich  has  been  used  to  fit 
scattering  data  for  the  system  uVNi(OOl)  (Gerlach  and  Hulpke, 


Potential  for  on-top  He/OC-Ni (001 ) system  over  several 
points  on  the  (y-z)  plane.  Starting  from  first  line  on 

3a/4,  a;  where  a is  the  lattice  spacing. 


V(Ra,Z)/eV 


Figure  5.5:  Potential  foi  bridge  He/OC-Ni(OOi > system  oyer  several 

the  right,  the  points  correspond  to  y=0,  a/4 
3a M , a;  vhere  a is  the  lattice  spacing 


l,  a/2, 


similar  equation  to  Eq.  (5.5) 


(5.6) 


The  term  V (R  ) is  identical  to  the  one  for  the  helium  potential  but 
Micha,  1963a).  The  term  VS(H)  contains  a repulsive  ond  an  attractive 


VR)  ■ Vrep<z>  * Vattr(x,y'z)  (5'7a) 

Vrep<z>  ‘ As  (5‘7b) 

Vattr(x,J,’z>  ■ B(e*Pl-26(z-z,)l  - 2 exp|-B(z-Zj)|)  <5. 7c) 

* 2cd>|-26(z-Zj)l  h(x,y)  (5.7c) 


5.1.  Hovever,  one  should  met 


is  already  being  explained  in  Sect 


term,  at  sufficiently  small  distances  z ft 
cancelled  by  the  repulsive  part  included  i 
parameters  used  for  the  Li*/0C-Ni(0010)  sy 


. The  potential 


again  found  to  be  steeply  repulsive  around  the  energies  of  interest 

potentials  for  the  top  and  bridge  model  respectively.  One  can 
observe  that  the  trends  for  turning  points  and  veil  depths,  vhen 
compared  between  models,  are  the  same  as  the  ones  found  for  the  He 

potential  in  the  y-z  plane  can  he  found  vhen  one  compares  the  two 


Table  5.2  Interaction  Potential  Parameters  for  Li*/OC-Ni(OI>l) 


1777.! 


As/CZ'BHjA 


a3/(z'°h)a 


Z/ Angstrc 


He/C0-Ni(001)  system, 
this  is  the  opposite  o: 


es  5.8  and  5.9  for  the  l.i*/0C-Ni <0M ) 
the  tunning  points,  as  one  moves  away 
hange  as  much  as  they  do  in  the 
is  because  the  M*-Ni(001)  interaction  i; 

occurs  Cor  the  He-Ni(OOl)  i 


ECfective  Classical  Difierential  C 


s-Sections 


The  interaction  potential  in  the  system  atom(ion)-ad: 

the  surface.  This  represents  a two  center  interaction  problem,  which 

classical  cross-sections  can  not  be  evaluated  from  the  formulation 
for  isotropic  potentials. 

For  anisotropic  potentials,  the  differential  classical  cross- 
sections  can  be  evaluated  from  a classical  trajectory  ?(t)  determined 


To  obtain  the  differential  cross-sections  one  needs  to  go  from  an 
initial  element  of  area  determined  by  the  impact  parameter,  b,  and 
the  initial  aaimuthal  angle,  to  a final  element  of  area 
determined  by  the  final  polar,  0,  and  aaimuthal.  *,  angles. 

Therefore,  one  needs  the  transformation  from  the  set  (b.d^l  to  the 
set  (0, 4) . This  transformation  is  obtained  by  the  Jacobian  and  the 
differential  cross-section  (Vilallonga  and  Micha,  1987)  is  then  given 
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Figure  5.8:  Potential  for  on-top  LiVOC-M(OOl)  system  over  several 

• points  on  the  (y-z)  plane.  Starting  from  solid  line  on 

.a;  vhere  a is  the  lattice  spacing. 


Figure  5.9:  Potential  for  bridge  Li*/0C-Ni(001)  system  over  several 

a;  where  a is  the  lattice  spacing. 


2j=  = (b/sinS)  ■ 


d inelastic 


molecular  differential  cross-sections  (Gentry,  1979). 

scattering  of  a beam  of  identical  particles,  with  a specified  kinetic 
energy,  Ej,  striking  a sample  of  tile  material  under  study.  The  atoms 

Using  a bundle  of  trajectories,  the  density  of  the  incident  beam 

represented  by  a grid  of  (x.,yj)  points.  These  points  are  related  to 
the  initial  conditions,  (b,&.),  of  the  trajectories  by  x . b cos$() 
and  y ® bsin$Q.  Denoting  by  dN^  the  number  of  trajectories  with 
final  scattering  angles  within  the  intervals  (6»d0)  and  (9±d$),  the 
differential  cross-sections  may  then  be  obtained  from 


do/di)  - (dNf/sin8  da  d»)/(dN./dxdy)  (9.10) 

The  quantity,  dM^,  can  be  extracted  from  studies  of  trajectories  by 
analyses  of  the  final  polar  and  azimuthal  angles  of  the  individual 
trajectories. 


Trajectories  vere  obtained  by  numerically  solving  Hamilton's 
equations  in  the  coordinate  system  shown  in  Pig.  5.1.  The  set  of  six 
first  order  differential  equations  was  integrated  by  a method  which 

extrapolation  (Shampine  and  Gordon,  1975).  The  initial  conditions 
(b,$Q,  Rmax)  were  chosen  by  requiring  that  the  interaction  potential 
be  small  with  respect  to  the  relative  translational  energy,  e.g.,  by 

The  area  of  the  target  sampled  by  the  trajectories  is  related  to 
the  maximum  impact  parameter,  bmax,  defined  as  the  largest  impact 

adsorbate.  The  area  was  chosen  to  be  given  by  flxfly,  vhere  Ox  = x|nax 
- Xq,  and  xfflax  a bBax  cos$q  (bmax>‘  Here  *o^raax^  *s  t*le  ' 11 ' ’ ‘ a ^ 

interaction  between  the  projectile  and  the  adsorbate  is  of 

three  separate  regions.  The  first  region,  for  small  impact 

strong  repulsive  interaction  with  the  adsorbed  molecule.  These 

second  region  is  where  the  trajectories  are  deflected  tvice.  The 
first  deflection  is  caused  by  the  interaction  with  the  adsorbate  and 

the  interaction  is  mainly  due  to  the  effect  of  the  surface.  These 


trajectories  have  impact  parameters 
in  the  calculation  and  vas  kept  to  • 

effect  of  the  image.  The  image  efft 
significance  for  large  deflections. 


tout  21!.  Several  of  the 

be  significantly  affected.  This  can  be 
trajectories,  in  the  (x-z)  plane,  for  the 


Li+  with  10  eV  of  energy  are  shovn. 

The  trajectories  show  the  shape  and  size  of  the  adsorbate  and,  as 
vas  expected,  they  reflect  the  different  sizes  of  the  adsorbate  in 
the  two  models,  see  figure  5.11. 

In  the  case  of  helium  atoms,  the  trajectories  end  at  larger  polar 


atoms,  and  trapping  does  not  occur.  Thir 
the  large  attractive  potential  in  the  cat 

The  turning  points  are  much  larger  ft 
lithium  ions.  Compare  Figures  5.10  and  5.12.  This  i 


f helium  it 


The  dependence  of  the  polar  scattering  angle  on  the  impact 
projectiles  and  both  models.  This  minimum  corresponds  to  the  impact 


Li*  trajectories  in  (x-z)  plane,  on-top  model;  for 


function  reaches  a maximum  which  indicates  the  beginning  of  the  third 
region  where  the  interaction  is  mainly  between  the  projectile  and  the 
surface.  These  features  are  shown  for  the  Li  ions  in  Figure  5.13 


The  differential  cross-sections  for  the  on-top  model  are 

of  these  results  is  that  the  scattering  is  predominantly 
backscattering.  In  the  case  of  Li  * ions,  double  collisions  represent 
the  main  contribution  for  all  scattering  polar  angles.  The  ratio  of 
single  to  double  contributions  depends  on  the  azimuthal  angle.  For 
lithium  ions,  the  single  collision  contribution  is  the  same  for 
different  directions  on  the  x-y  plane,  but  the  double  contribution  is 
not  symmetric  with  respect  to  the  z axis.  In  the  case  of  helium 
atoms,  both  contributions  are  angle  dependent.  A possible 
explanation  for  this  is  the  larger  corrugation  present  in  the  helium 
case.  Also,  one  can  notice  that  at  certain  azimuthal  angles,  e.g. 
♦=45°,  the  importance  of  single  collisions  with  respect  to  double 
ones  is  reversed  for  certain  ranges  of  polar  angles.  This  factor 
arises  from  the  symmetry  of  the  model.  That  is,  the  nearest  neighbor 
distance  along  the  y axis,  *-90“,  is  approximately  2.489  A,  whereas, 
in  the  direction  with  *-45°,  it  is  3.52  A,  thus  giving  a larger 


- 0”( — ), 


b/Angstroms 


He  deflection  function,  on-top  model.  = 0#< ), 

45°( — ).  The  inset  shoes  an  enlargement  Eor  large 
impact  parameters. 


Figure  5.16:  He  differential  cross-section,  do/da,  on-top  model,  for 
single  ( ) and  double  ( — ) collisions;  a)  6 - 0B  b) 


CHAPTER  6 

SY  TRANSFER  IN  HYPERTHERMAL  COLLISIONS 
OP  He  AND  Li'  WITH  CO  ADSORBED  ON  Ni(OOl) 

In  this  chapter  we  apply  the  theoretical  model  presented  in 
Chapter  A to  the  scattering  of  He  and  Li*  from  a CO  molecule  adsoi 
on  a Ni(OOl)  surface.  Ve  calculate  double  differential  cross- 
sections  for  the  energy  and  angular  distributions  of  the  He  atoms 
the  Li  ions.  Results  are  obtained  for  experimental  conditions  ol 
high  relative  energy,  1 eV  for  He  atoms  and  10  eV  for  Li',  for  noi 
incidence  and  for  the  low  coverages  where  the  theoretical  model  i 1 
applicable. 


6.1  Double  Differential  Cross-Section  for  the  System  He/OC-Ni(001 
The  calculations  are  based  on  the  expressions  obtained  in  Chapter 
A using  the  short  time  expansion  of  the  vibrational  TCFs  for  the  low 
frequency  modes  of  the  target  which,  in  our  case,  are  all  the  surface 
and  cluster  modes  except  the  CO  stretch. 

According  to  Eq.  (A, 22),  the  partial  double  differential  cross- 

will  be  given  by  a Gaussian  distribution  of  energy  transfer,  with  a 

for  both  mechanisms  we  are  considering  here;  that  is,  for  single  and 

We  have  done  the  calculations  considering  single  and  double 
collisions.  We  find  that  single  collisions  arc  the  predominant 
mechanism  for  helium  scattering.  In  figures  6.1  and  6.2,  the  partial 


the  single  collision  contribution  is  larger  than  the 


and  -n-2 , respectively,  vhere  it  also  can  be  seen  that  single 


collisions  are  even  more 
transferred  into  the  CO 
Double  differential 


cross-sections  (DDCS)  at  scattering  angles  of 
, for  the  top  and  bridge  models  in  Figures  6.5 


increasing  e,  vhich  means  vith  increasing  E , since  E,  is  fixed  by  n. 
As  E is  further  increased,  the  DDCS  goes  through  a maximum  and  then 
decreases  to  zero.  In  the  case  of  no  excitation  of  the  CO  stretch. 


that  is,  for  -n»0,  there  is  no  apparent  difference  betveen  the  tvo 

bridge  model  gives  larger  quantities  for  the  DDCS  and  Gaussians 
centered  at  slightly  smaller  values  of  transferred  energy.  This  is 
found  to  be  in  accordance  vith  the  fact  that  the  frequency  of  the  CO 


stretch  in  the  bridge  model  is  slightly  smaller  than  in  the  top  case 
Gaussian  is  given  by  (£  * A) . 

increases  when  the  polar  scattering  angle  increases.  This  trend  is 
easily  understood  vhen  one  recalls  that  the  CO  stretch  is  a mode  vith 


Aj  symmetry  and  motions  perpendicular  to  the  surface.  Therefore, 
collisions  at  small  impact  parameters,  leading  to  backyard 
deflections,  vould  be  the  most  efficient  to  excite  the  CO  stretch. 

It  is  also  seen  that  the  same  trend  is  followed  when  n=0 . In  this 
case  the  overall  process  involves  the  simultaneous  excitation  of  the 
slow  modes,  some  of  vhich  are  E type  (motions  parallel  to  the 


and  some  of  which  are  A type  and  more  likely  to  be  exci 
trends  compound  to  give  the  resulting  angle  dependence 
indicates  that  the  A type  modes  are  the  dominant  ones. 
All  the  results  presented  and  discussed  up  to  nou  u 


obtained 


rmined  by  its  magnitude 
e frequencies  satisfy  ho 


s expected,  the  peaks 


ition  to  the  quantity  K_T.  Modes 
will  have  excitation  probabilities 
ik  , Independent  of  temperature  for 

whose  frequencies  satisfy  Mw-S  k„T,  and  is 
of  the  Gaussians-  From  equation  (4.22)  one 

Figures  6.10,  6.11,  and  6.12,  we 

•e  narrower  than  the  ones  presented  in 


calculation 


snail  peak  appears  in  the  region  of  very  small  energy  transfer.  In 
Figure  6.13  this  peak  is  isolated  and  found  to  be  centered  at  a point 
vhere  the  energy  transferred,  e is  about  17  meV,  and  presents  a 


Poisson,  rather  than  a Gaussian,  distribution.  This  seems  to 
indicate  that  it  corresponds  to  a mode  with  a frequency  satisfying 
the  minimum  condition  of  M«.  > 8 meV.  In  our  analyses  of  normal 

frequencies  of  15. 5 meV  and  16.2  MeV  in  the  bridge  and  top  cases, 
respectively.  This  peak  may  then  be  associated  vith  these  modes, 
especially  considering  that  it  appears  in  both  models. 

Finally,  the  dependence  of  the  DDC5  vith  the  azimuthal  angle  6 

The  trends  followed  by  the  two  models  are  essentially  the  same,  with 
a smaller  DDCS  in  the  (110)  direction  (9=0,  90*)  than  in  the  (100) 
direction  (0-45°).  This  result  indicates  that  the  projectile  does 
not  seem  to  notice  the  corrugation  of  the  surface. 


6.2  Double  Differential  Cross-Sections  for  the  System  Li*/0C-Ni(001) 
In  the  calculations  for  this  system,  we  have  chosen  energies  in 
the  range  of  10  to  20  eV,  which  is  considered  an  intermediate  regime 
between  low  energy  (thermal)  and  high  energy,  or  Rutherford 
scattering,  vhere  the  dynamics  are  rather  simple.  The  calculations 
are  based  on  the  same  formulation  used  for  the  He/OC-Ni(001).  As 
before,  we  consider  the  separation  of  the  modes  into  fast  ones,  that 


CO  stretch, 


more  energy  is  available  to  be  transf 
the  fast  mode  must  be  considered, 
as  a function  of  energy  transferred  f 


e differential  cross-sections 


and  4=0°  in  Figure  6.16. 
for  He  atoms,  two  peaks, 
peak  at  high  energies  is 
mechanism  and  the  one  al 


I to  originate  from  the  double  collision 

As  more  energy  is  transferred  into  the  fast  mode,  that  is,  for 
increasing  values  of  n,  the  single  collision  starts  to  become  more 
important.  At  this  value  of  9,  the  excitation  of  the  (X)  stretch  by 

always  the  case  can  be  seen  from  Figures  6.17,  6.18,  and  6.19,  in 
which  similar  results  are  represented  for  9-110,  150,  and  170®, 
respectively.  As  one  goes  from  small  to  large  6,  the  probability  of 
exciting  the  CO  stretch  by  a double  collision  mechanism  becomes 
smaller  at  the  same  time  that  the  probability  of  exciting  the  slov 
modes,  represented  by  n=0,  becomes  larger  by  the  same  mechanism. 

This  trend  could  be  explained  if  one  tries  to  visualize  the 
trajectory  of  a double  collision  with  a large  scattering  angle.  Such 
a trajectory  would  be  one  in  which 
small  deflection  by  the  adsorbate  e 
surface  with  a very  small  angle  of 


e projectile  In 


Consequently, 
projectile  vould  be  deflected 


modes  of  the  metal  surface.  The  fact  that  the  single  collision  peak 
almost  disappears  in  going  from  8=100  to  0=170°  seems  to  suggest  that 

excitation  than  the  modes  with  A symmetry,  giving  lower  DDCS  for 
larger  scattering  angles.  The  total  DDCS  obtained  by  summing  over 

trend  observed  in  the  individual  contributions  is  seen  here. 

Results  obtained  for  the  bridge  model  are  not  much  different  than 
the  ones  obtained  for  the  top,  although  some  differences  can  be 


Additionally, 


lmost  the  same  probabilities, 
e bridge  model  appear  to  be  vie 


larger  number  of  modes  considered  as  slow  modes  in  the  bridge  model. 
In  Figure  6.22  a comparison  of  the  DDCS  summed  over  all  n for  the  two 
models  is  presented  for  ©-100°.  It  is  found  that  the  bridge  model 
gives  smaller  DDCS  than  the  top  model.  The  effect  of  the  temperature 
is  shown  in  Figure  6.23  where  the  DDCS  summed  over  all  n are 


presented  for  a temperature  of  100°  K.  As  in  the  case  of  scattering 


by  He  atoms,  the  width  of  the  peaks  decreases  when  decreasing  the 


collision  peak 
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SUMMARY 


CONCLUSIONS 


In  this  work  we  have  presented  a theoretical  model  suitable  for 
the  study  of  scattering  of  atoms  or  ions  with  hyperthermal  energies 
by  molecules  adsorbed  on  metal  surfaces  at  low  coverages.  A simple 
model  to  obtain  the  vibrational  modes  of  molecules  adsorbed  on 
surfaces  has  been  presented  and  frequencies  have  been  calculated  for 
a CO  molecule  adsorbed  on  a Ni(001)  surface.  The  procedure  allows 
for  the  modelling  of  different  adsorption  sires  and  makes  the 
classification  of  the  modes  according  to  symmetry  an  easy  task. 

An  interaction  potential  has  been  constructed  and  a trajectory 
study  of  the  elastic  scattering  of  He  atoms  and  Li*  ions  by  a CO 
molecule  adsorbed  on  a metal  surface  has  been  carried  out.  It  is 


10  eV  for  Li  atoms,  the  effective  classical 
adsorbed  molecule-surface  complex  do  not  shot 

potential  effect  is  almost  negligible  except 


Ve  have  obtained  double  differential  cross-sections  as  functions 
of  energy  transfer,  scattering  angles,  incident  energy  and  metal 
temperature  for  the  cases  of  a CO  molecule  adsorbed  either  at  the 
bridge  or  on-top  configuration. 


excitation  process  is  the  most  likely  to  occur  at  a given  scattering 

predict  that  in  the  scattering  of  He  atoms,  the  most  likely  energy 
loss  vould  be  found  at  0.4  eV  and  vould  be  due  to  a collective 
excitation  of  modes  vith  lover  frequencies  and  not  to  the  excitation 


large  transfer 
excitations  of 


iso  say  that  features  found  at 
? adsorbed  molecule. 


the  excitation  of  individual  modes  can  be  obtained  by  performing 
experiments  at  lover  temperatures  where  features  of  the  lov  frequency 


low  CO  coverages  on  a Pt(001)  surface,  i 
attributed  t. 


molecules  (lahee,  i 


:0  vibrations.  Another  very 
16.5  meV  and  was  assigned  to  the  hindered  rotation.  This  assignment 
was  based  on  the  results  of  some  calculations  (Richardson  and 
Bradshaw)  that  have  given  a value  of  22. 8 meV  for  this  mode. 

Hovever,  in  our  calculations  of  the  vibrational  modes  we  obtained 
values  for  the  frustrated  rotations  for  the  bridge  and  top  models,  of 

centered  at  about  17  meV  for  both  models.  Ue  believe  this  peak 
should  correspond  to  the  frustrated  rotations  on  either  model. 


For  one  thing,  double  collisions  seen 
process.  Although  excitations  of  the 

predictions  can  be  made.  For  example, 


i ions  suggest  t 
comparison  vith  the  He  case, 
i to  play  a major  role  in  the 
) CO  stretch  involve  large 
they  give  results  vhete 

found  at  large  energy  transfers, 
the  scattering  angle,  one  can  exp. 


another  narrov  peak  centered  at  larger  va.' 

more  detailed  calculation,  where  more  modes  are  specified  by  the 
quantum  numbers  n,,  is  needed.  This  can  be  readily  obtained  from  01 
approach  when  the  short  time  expansion  approximation  is  only  used  ft 
a few  modes.  For  example,  the  bridge  site,  where  the  distinction 
between  high  and  low  frequency  modes  is  not  so  sharp,  could  be 


A more  rigorous  treatment  of  double  and  multiple  collisions  can 
be  based  on  the  multiple  scattering  expansion  of  the  transition 

One  obvious  limitation  in  our  model  is  that  it  ignores  the  large 

These  modes  could  in  principle  be  obtained  with  the  standard  methods 
of  lattice  dynamics,  and  have  actually  been  used  in  the  problem  of 
adsorbate  vibrations  (Rahman  et  al.,  1982).  However,  these  methods 


applied  to  molecular  adsorbates  vi i 


since  they  rely  on  lattice  periodicity  their  application  to 
disordered  systems  is  even  more  complex. 

A possible  method  which  could  be  used  to  deal  uith  this 
limitation,  but  vould  still  have  the  convenience  of  a cluster  model, 
consists  of  including  periodic  boundary  conditions  for  the  clusters. 
Periodicity  in  the  directions  parallel  to  the  crystal  surface  is 
accomplished  by  adding  new  bond-stretching  internal  coordinates,  in 


additional  to  the  ones  within  the  cluster,  betweer 

two  adjacent  layers  of  the  crystal.  This  method  t 

allowing  thus  to  examine  the  dispersion  of  the  moc 
energy  transfer  into  these  modes  as  a function  of 

O/Ni {001 ) system  (Lloyd  and  Hemminger,  1985)  and  f 
almost  exactly  the  results  obtained  with  the  Greer 


edge  atoms  located 
:n's  function  lattice 


hyperthermal  energies  which  are  dominated  b 
these  modes  are  not  likely  to  have  an  effect  in  the  process. 

assumption  of  harmonicity.  One  of  the  most  interesting  aspects 
associated  uith  the  vibrational  states  of  chemisorbed  molecules  i: 
the  change  in  vibrational  frequencies  and  lifetime  of  a molecule  i 

broadening  of  the  vibrational  state  of  a chemisorbed  molecule  on  . 
metal  surface  is  vibrational  energy  relaxation.  This  mechanism  i: 
associated  with  the  coupling  between  the  localized  vibrational  mo* 
of  the  adsorbated  and  the  phonon  modes  of  the  metal.  The 


potential  and  also  on  the  frequency  of  the  localized  vibrational  mode 
in  relation  to  the  highest  phonon  frequency  of  the  metal.  By  analogy 


1980)  one  may  infer  that  for  modes  with  frequencies  greater  than 
twice  the  highest  phonon  frequency  of  the  metal  the  damping  provided 


lov-frequency  modes,  such  as  the  C-Ni  stretching  vibration  fi 
adsorbed  on  a Ni(001)  surface,  decay  via  tvo-phonon  emission 
energetically  possible  and  probably  the  relevant  energy  rt 
process.  Again  because  of  the  short-time  dynamics  in  hyperthermal 

If  necessary,  at  low  collision  energies,  anharmonicity  can  be 
introduced  in  the  formalism  of  the  time-correlation  functions  for 

the  case  of  polyatomics  (Vilallonga  and  Micha,  1983b). 

In  a general  way  one  can  write  the  vibrational  potential  in  term: 

of  the  normal  coordinates  that  diagonalize  the  harmonic  part  of  the 
vibrational  Hamiltonian,  the  total  vibrational  Hamiltonian  is  given 


Hh  ■ <ETE  * 9TSh  0)/2  , 

vhere  P is  the  column  matrix  of  momenta  P.  - -iha/3i). . One  then 
finds  that  the  TCFs  are  given  in  term  of  displacement-displacement 
correlation  functions  (DDCS),  «0(0)  Q1(t)»,  which  can  be  evaluated 
from  the  corresponding  double-time  Green  functions.  The  procedure 


" - *,D  . S*  vi  ih  DTD  - 1 and 


Hg  cothj  . 2DTU  D . 

W - «Q  QT»  ^ - «0  £T»  , 


V*-'*  ■ n “IX-Xj  "US,  * 


* X.Iexpf-ij  - ly)  . exp (ij  . iij t)|/2}  , 
>j  ■ bll-ljAai,  slnhij)  , 
c-CD. 


H(t), 


|«,«>  - expl-iHjCt-.jVKI  Wl.tj)  |«.ti> 

"i  . li.  ii(t) 


ar-  2 **.«- 


dP/de  = J dt/2nh  F(t)  exp(-ict/H) 


Ux(t)  = exp(iHxt/H)  expfii 


the  integration  o 


energy  transfer  is  given  by 

■ semiclassical  evolution 

■ will  involve:  first  defining 


H'(titi)  = 0K(t)  . exp| i(Hi(t-ti)/H|  V(t)  exp|-iH.(t-t.)/H| 


dK(t)  - |?2<t)  - ?jl/(2M) 

f being  the  initial  relative  momentum.  This  approach  would  provide 
the  means  to  tackle  the  problem  of  scattering  of  atoms  at  thermal 


energies  by  adsorbates. 


collisional  TCFs 


general  validity. 


special  cases.  Several  aspects  of  the  approach  that  ve  have 
presented  can  be  developed  in  future  vork.  In  particular,in  view  of 
the  increasing  amount  of  experimental  work  currently  in  progress  on 
the  scattering  of  atoms  at  thermal  energies,  it  vould  be  fruitful  to 


look  into  the  semiclassical  method. 

applied  without  major  changes  to  different  systems.  Calculations 
could  be  done  for  example  for  a different  crystal  face  of  the  same 
metal  or  any  other  metal,  or  for  a different  type  of  adsorbed 
molecule  and  even  for  a different  type  of  projectile.  Also, 
calculations  could  be  done  for  incident  angles  different  from  zero. 
The  case  of  higher  converages  could  also  be  studied  with  the  present 
approach,  including  the  interaction  between  the  adsorbates. 
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'<JtNi  * YN1)/2  * **  ZNi/2 

-<Xm  - Yn,)/2  * 42  Zm/2) 
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ridge  Cluster 
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Appendix  A (continued) 
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